Alternating hexagonal and striped patterns in Faraday surface waves 
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A direct numerical simulation of Faraday waves is carried out in a minimal hexagonal domain. 
Over long times, we observe the alternation of patterns we call quasi-hexagons and beaded stripes. 
The symmetries and spatial Fourier spectra of these patterns are analyzed. 



The Faraday instability [l[ describes the generation of 
surface waves between two superposed fluid layers sub- 
jected to periodic vertical vibration. Although these 
waves usually form crystalline patterns, i.e. stripes, 
squares, or hexagons, they can form more complicated 
structures such as quasicrystals or superlattices [2|-|5[ . We 
have recently carried out the first three-dimensional non- 
linear simulations of the Faraday instability Q and have 
reproduced the square and hexagonal patterns seen in 
[?| with the same physical parameters; see Fig. [TJ How- 
ever, the hexagonal pattern is not sustained indefinitely, 
but is succeeded by recurrent alternation between quasi- 
hexagonal and beaded striped patterns. This long-time 
behavior is the subject of the present study. 

The first detailed spatio-temporal experimental mea- 
surements of the interface height z = C,(x,y,t) were un- 
dertaken by Their optical technique required the 
two fluid layers to have the same refractive index, which 
led them to use fluids of similar viscosities and den- 
sities: pi = 1346 kgm" 3 , i/i = 5.35 x 10 -6 m 2 s -1 , 
p-2 = 949 kgm -3 , V2 = 2.11 x 10~ 5 m 2 s _1 and surface 
tension a = 35mNm _1 . These parameters, especially 
the density ratio pil P\ = 0-7, differ markedly from most 
studies of Faraday waves, which use air above either wa- 
ter or silicone oil and so have pij p\ ~ 0.001. At rest, 
the heavy and light fluids occupy heights of h\ = 1.6 mm 
and /12 = 8.4 mm, respectively. The imposed vibration 
has frequency / = 12 Hz and the Faraday instability 
leads to subharmonic standing waves, so that Q{x, y, t) 
oscillates with period T = 2/f = 0.1666 s. Floquet 
analysis Q for these parameters yields a critical wave- 
length of A c = 2ir/k c = 13.2 mm, with which the exper- 
iments show close agreement 0- Thus, another atypical 
feature of this parameter regime is that h\ <C A c ; see 
Fig. [TJ Floquet analysis also yields a critical acceleration 
of a c = 25.8 ms~ 2 = 2.63 g. Our simulations are carried 
out at a = 38.0 ms -2 = 3.875 g = 1.473 a c , for which 
hexagons were observed experimentally Q- 

We summarize our formulation and the numerical 
methods used to compute the fluid motion; see @ for 
a more detailed description. Our computations use a 



FIG. 1. Above: instantaneous realization of hexagaonal Fara- 
day waves in a domain of size L x x L y x (hi + h?)- The 
density and viscosity are p2, f2 above the free surface at 
z = ((x,y,t); below, they are pi and v\. The displacement 
amplitude of ( is comparable to hi; at its lowest points, ( 
almost touches the lower boundary. Below: projection of the 
horizontal doubly periodic rectangular domain of dimensions 
L x x L y — 2A c /v3 x 2A C shows that it is compatible with 
a hexagonal lattice. The minima of £ (above) are located at 
the vertices and centers of the hexagons (below). The right 
corner represents a portion of the rectangular computational 
grid; the actual grid is twice as fine in each direction. 



single-fluid model, representing the velocity u and pres- 
sure p over the whole domain on a staggered MAC mesh 
Q which is fixed and uniform. The viscosity and den- 
sity are variable, taking the values v\, p\ for the denser 
lower fluid, vi, pi for the lighter upper fluid and vary- 
ing over a few gridpoints at the interface. The mov- 
ing interface, defined by z = £(x,y,t), is computed by 
a front-tracking ficl] / immersed-boundarv 11 1 method on 
a semi-Lagrangian triangular mesh which is fixed in the 
horizontal x and y directions and moves along the verti- 
cal direction z. The interface is advected and the density 
and viscosity fields updated. The capillary force is com- 
puted locally on the Lagrangian mesh and included in 
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FIG. 2. Maximum interface heights max^y (,(x, y, t) (rapidly oscillating curve) and mayi x ^y^ t ^+T}C,{x,y,t) (smooth envelope). 
Surrounding visualizations show instantaneous contour plots of C,{x,y,t) at times indicated by black dots on rapidly oscillating 
curve and along abscissa. Red/blue colors regions of maximum/minimum interface height. The color map of each contour 
plot is scaled independently of the others. Over the large white areas, the interface is very close to the bottom and almost 
flat. Visualizations shown at ti + jT/4 for j — 0, 1,2,3 (i.e. over one subharmonic period) for hexagons at ti, symmetric and 
nonsymmetric beaded stripes at ti and £4 and quasi-hexagons at tz and £5. 



the Navier-Stokes equations, which are solved by a pro- 
jection method. The computations are carried out in 
the oscillating reference frame of the container by adding 
a time-periodic vertical acceleration asin(27r/i)e z to the 
equations of motion. No-slip boundary conditions are im- 
posed at the top and bottom boundaries, while periodic 
boundary conditions are used at the vertical boundaries. 

The horizontal dimensions of the domain are chosen to 
accomodate a hexagonal pattern. We take L x = 2A c /-\/3 
and L y = 2A C , as shown in Fig. [TJ so that large-scale 
spatial variations are inaccessible. This domain is also 
compatible with striped or rectangular patterns, as will 
be discussed below. The simulations were run with a spa- 
tial resolution of N x x N y x N z = 58 x 100 x 180. Each 
horizontal rectangle is subdivided into 64 triangles to rep- 
resent the interface. To validate the spatial discretiza- 
tion, we repeated the simulations with a finer resolution 
of N x x N y x N z = 75 x 125 x 225. Although small quanti- 
tative changes were seen, the dynamics remained qualita- 
tively unchanged. The timestep is limited by the advec- 
tive step, taking values varying between T/24000 for a 
hexagonal pattern and T/4000 for a beaded striped pat- 
tern. This makes the simulation of behavior over many 
subharmonic periods extremely time-consuming. 

Starting from zero velocity and an initial random per- 
turbation of the flat interface, our simulations produced a 
hexagonal pattern which oscillates subharmonically with 
the same spatio-temporal spectrum as , as detailed in 
Q. In the experiments, hexagons are transient and dif- 



ficult to stabilize 1121 . competing with squares and dis- 
ordered states fH Il3|. In our simulations, after about 
10 subharmonic periods, we observed a drastic departure 
from hexagonal symmetry. Figure [2] shows the instanta- 
neous maximum height max^ . y £(x, y, t) and its envelope 
max x t y t [t t t + T]^{x,y,t). Surrounding the time-evolution 
plot are contours of the interface height at representative 
times over one subharmonic cycle. 

At times t\ + jT/4, the patterns are hexagonal. Each 
is invariant under rotations by 7r/3 and under reflection, 
for example about y = nX c ; these operations generate 
their isotropy subgroup (group of symmetries), which is 
isomorphic to Dq. We call the patterns at ti + jT/4 
beaded stripes. They satisfy the two symmetry relations: 

C(x, n\ c -y) = CO, y) = C(mXc/V3+x -x, y+n\ c ) (1) 

where xq sa A c /(2v / 3) is a spatial phase. The second 
equality in ([T]) is variously called shift-and-rcflcct or glide- 
reflection symmetry. These invariances describe the crys- 
tallographic group called pmg or p2mg [13l . |14| , which is 
isomorphic to Z 2 x Z 2 . At later times, the patterns have 
no exact symmetries. Nevertheless, the patterns at times 
£3 + T/4 and £5 + 3T/4 contain large flat cells surrounded 
by six small peaks like those at t\ + 3T/4, which lead us 
to call them quasi-hexagons. Quasi-hexagons at £3 and 
£5 appear in two distinct forms, which are related by the 
spatio-temporal symmetry 

C(mA c /V3 + xo - x, y + y , t 3 + T/2) = ((x, y, t 5 ) (2) 
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FIG. 3. Time-filtered spatial Fourier spectra £mn(ti): hexagons at ti, beaded stripes at t2, quasi-hexagons at £3 and t$, 
nonsymmetric beaded stripes at t± and t&. 



with phases xo ~ 0.7X c /^/3 and yo w A c /2. The patterns 
at *4 and ^ resemble those at ti but do not satisfy ([TJ; 
we call these nonsymmetric beaded stripes. These are 
also related by the spatio-temporal symmetry ©. 

To quantify the behavior in Fig. [5J we have studied 
the spatial Fourier spectrum. The rectangular box of 
Fig. [1] constrains the wavevectors to lie on the grid k„ m = 
[\fime, x + ne y ) k c /2, m 7 n <G Z, as shown in Fig. [3] and 
H^,. We define the time-filtered spatial Fourier transform: 

C(x,i) 



A, 



X ( mn (t) , (mn(t) = ^ax^ |Cmn(*)| 



In Fig. [31 we plot Cmn(*i) for ti . . .t 6 . For the hexago- 
nal pattern at time t\ the modes with non-negligeable 
amplitude are those belonging to hexagons, primarily 
(m,n) = (±1,±1) and (0,±2), with k/k c = 1, but also 
(±2,0) and (±1,±3), with k/k c = and (0,±4) and 
(±2, ±2), with k/k c = 2. The spectrum at ti is dom- 
inated by modes (0, ±2) and (±1,0), which combine to 
form the beaded striped patterns - with one bead over L x 
and two stripes over L y - seen in Fig. [5] at ti. The spec- 
tra at tj, and t§ combine wavenumbers seen at t\ and ti, 
while those at £4 and t§ are asymmetric versions of that 
at i 2 - The spatio-temporal symmetry ^ is manifested 

byC (^3,4) — Cmn{t5fi)- 

Figure 0b shows the time evolution of ( mn (t) for the 
dominant wavevectors, grouping those with the same k 
value into a single graph. Because Cm.-n = C—m.n (result- 
ing from the reality condition ( m -„ = (l mn ), we plot 
only positive n. We first describe the evolution of the 
four most dominant modes in the leftmost column. Dur- 
ing the hexagonal phase at t±, the only modes with non- 
negligcablc amplitude are those belonging to hexagons. 
At t/T sa 40, these drop abruptly, with the exception of 
(02 ■ This is accompanied by a burst in amplitude of some 
of the non- hexagonal modes, notably (jo, followed by its 
saturation. Shortly after ti, Ci,i begins to rise, followed 
by £-1,1) eventually C10 rises as well; we conjecture that 
its growth is fueled by mode interactions arising from 
kn ±k_i,_i + ki,o = ki . By t 3 , these have attained val- 
ues which approach £02, leading to the quasi- hexagonal 
patterns at £3. Mode amplitudes C±i.i then fall quickly, 
followed by £1.0, leading again to a short-lived beaded 



striped pattern at time £4. The cycle then repeats, but 
this time C-1.1 rises before £1.1 and attains a higher peak 
at time £5, leading to the difference between the quasi- 
hexagonal patterns at £5 and £3. The next cycle shows di 
leading again at t&. Most of the higher k modes behave 
like their lower k analogs. 

Figure 3J; shows a phase portrait, projecting the dy- 
namics onto coordinates, (11 + C-i,ij C11 — C-1,1 an d 
£10 to represent the dynamics. Here, the simulation has 
been extended and shows several additional repetitions. 
The concentration of points indicate that the hexago- 
nal pattern at time t\ and the beaded striped pattern 
at t<x arc saddles. Afterwards, the trajectory consists of 
two crossed loops connecting £3, £4, £5 and t§. Several 
dynamical-systems scenarios lead to limit cycles which 
visit symmetrically related sets, e.g. Hopf bifurcations 
[TBI ] or heteroclinic cycles Since the saddles at t\ 

and ti are not part of the periodic cycle, nor its center, 
our planned investigation of its origin will be challenging. 

A complete analysis of the general bifurcation problem 
on a hexagonal lattice 17 , HI shows that the patterns at 
onset are hexagons and stripes. This analysis of steady 
states on a hexagonal grid does not apply directly to our 
study of subharmonic standing waves on the rectangular 
lattice of Fig. [TJ At the other extreme, in large-domain 
experiments, hexagonal patterns are often observed to 
undergo other spatio-temporal dynamics, such as com- 
petition with squares @, 0, EH , which are inaccessible 
to current mathematical analysis and to our simulation. 

We have observed complex long-time temporal behav- 
ior in a fully resolved three-dimensional simulation of 
Faraday waves in the minimal domain which can acco- 
modate a hexagonal pattern. Although this scenario may 
prove to be replaced by other dynamics in large domains, 
we believe that it is of interest in its own right and that it 
may well be applicable to other pattern-forming systems. 

We thank P.-L. Buono and M. Golubitsky for sharing 
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the Natural Sciences and Engineering Research Council 
of Canada. 
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FIG. 4. a) Spatial Fourier grid of the domain. Hexagonal symbols: sets of six k m „ which share the same length k. This portion 
of the grid contains three hexagons, with k/k c = 1, v3 and 2. The hexagonal modes are located at the vertices (and not the 
sides) of the large hexagons indicated by solid and dashed lines. Square symbols: sets of four k mn which share the same k. 
Circles: other modes. Colors differentiate between k mn and its images under rotations of Itt/3 for hexagonal modes or under 
reflections across the n axis for the others, b) Temporal evolution of C, m n(t), grouped by length k. The color code is that 
used in the spatial Fourier grid. Dashed blue curves: hexagonal modes (0,2), (1,3), (0,4) and non-hexagonal modes (—1,2), 
(—2,1). Solid black curves: hexagonal modes (1,1), (2,0), (2,2) and non-hexagonal modes (1,0), (1,2), (2,1). Dash-dotted 
orange curves: hexagonal modes (—1, 1), (—1, 3), (—2, 2). c) 3D phase portrait of temporal evolution of modes for an extended 
simulation. Points are equally spaced in time; their color evolves in time from dark to light. Times t\ . . . t$ marked by crosses 
and by vertical lines in b). 
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